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Abstract 

^ ' Starting from kicked equations of motion with derivatives of non-integer orders, we 

(^ '. 

If^ , obtain " fractional" discrete maps. These maps are generalizations of well-known universal, 

0^■ 

ff^ ■ standard, dissipative, kicked damped rotator maps. The main property of the suggested 



o 



% 



fractional maps is a long-term memory. The memory effects in the fractional discrete maps 
mean that their present state evolution depends on all past states with special forms of 
weights. These forms are represented by combinations of power-law functions. 

1 Introduction 

There are a number of distinct areas of physics where the basic problems can be reduced to 
the study of simple discrete maps. Discrete maps have been used for the study of evolution 
problems, possibly as a substitute of differential equations [H |2l [3l S]. They lead to a much 
simpler formalism, which is particularly useful in simulations. The standard map is one of the 
most widely studied maps. In this paper, we consider fractional generalizations of discrete maps 
that can be used to study the evolution described by fractional differential equations [3 |6l [7]. 



The treatment of nonlinear dynamics in terms of discrete maps is a very important step in 
understanding the quahtative behavior of continuous systems described by differential equa- 
tions. Note that the continuous limit of discrete systems with long-range interactions gives 
differential equations with derivatives of non-integer orders (see for example, [El |9l [in])- The 
derivatives of non-integer orders are a natural generalization of the ordinary differentiation 
of integer order. Fractional differentiation with respect to time is characterized by long-term 
memory effects that correspond to intrinsic dissipative processes in the physical systems. The 
application of memory effects to discrete maps means that their present state evolution depends 
on all past states [m [121 HSl [12 [IS [16] . 

The mapping Xn+i = f{xn) does not have any memory, as the value Xn+i only depends 
on Xn- The introduction of memory means that the discrete value Xn+i is connected with the 
previous values x„, x„_i, . . . , Xi. Particularly, any system, which is described by a discrete 
map, will have a full memory [16], if each state of the system is a simple sum of all previous 
states: 

n 

Xn+l = ^f{Xk), (1) 

fc=l 

where f{x) is a function that defines the discrete map. In general, the expression of equation 
([1]) can tend to infinity. Note that the full memory exists for functions that give a finite sum 
in Eq. ([T]). The full memory is ideal because it has the same action upon the next states as all 
the others in memory. The map with a long-term memory can be expressed as 

n 

Xn+i = ^ Vain, k)f{Xk), (2) 

fe=l 

where the weights V^(n, fc), and the parameter a characterize the non-ideal memory effects. 
The forms of the functions f{x) and \4(n, fc) in equation ([2]) are obtained by the differential 
equation. Therefore the conditions on f{x) are not discussed. Note that linear differential 
(kicked) equations give the linear function f{x) = x. The discrete maps with memory are 
considered, for example, in the papers |Tll [121 [ISl [H [13 [16]. The important question is a 
connection of fractional equation of motion and the discrete maps with memory. It is important 
to derive discrete maps with memory from equation of motion with fractional derivatives. 



Here, we point out the well-known notions such that the full memory and the long-term 
memory (LTM). For the one-dimensional case, the full memory is defined by equation ([1]), and 
the LTM is defined by Eq. ([2]). In general, the expression of Eq. ([1]) can tend to infinity. Note 
that the full memory exists for functions that give a finite sum. This is a basic condition on 
f{x) in Eq. ([1]). In the paper, we consider the LTM for fractional equations. For the one- 
dimensional case, the LTM is defined by equation ([2]). The one-dimensional case is considered 
for simplification. In this paper, we consider a two-dimensional case of variables (x„,p„). Note 
that the conditions on the function f{x) of the LTM are not defined a priori. This function is 
obtained from an equation of motion. Therefore the conditions on f{x) are not discussed. The 
properties of the function /(x) are defined by the kicked differential equation. For example, 
the functions /(x), which are derived from linear differential (kicked) equations, are linear 
(fix) = x). 

It was shown [17] that perturbed by a periodic force, the nonlinear system with a fractional 
derivative exhibits a new type of chaotic motion called the fractional chaotic attractor. The 
fractional discrete maps allow us to study new type of attractors that are called pseudochaotic 

m- 

In section 2, a brief review of discrete maps is considered to fix notation and provide conve- 
nient references. In section 3, a fractional generalization of the universal map is obtained from 
kicked fractional equations of motion with order 1 < a < 2. We prove that the usual universal 
map is a special case of the fractional universal map. Some examples of fractional universal 
map are suggested. In section 4, a fractional universal map for the case of a > 2 is obtained. 
In section 5, a fractional kicked damped rotator map is derived. Finally, a short conclusion is 
given in section 6. 

2 Universal and standard map 

In this section, a brief review of well-known discrete maps is considered to fix notations and 
provide convenient references. For details, see [U El |3l H] . 



Let us consider the equations of motion 

oo 

x + KG{x)Y,s{f-n)=0 (3) 

n=0 

in which perturbation is a periodic sequence of delta-function-type pulses (kicks) following with 
period T = 27r/i/, K is an amplitude of the pulses, and G{x) is a some function. This equation 
can be presented in the Hamiltonian form 

x = p, p + KG{x)J2^{f-^)=0- (4) 

n=0 

It is well-known that these equations can be represented in the form of discrete map (see 
for example |4J). Between any two kicks there is a free motion 

p = const, X = pt + const. (5) 

The solution of the left side of the n-th kick 

Xn = xitn — 0) = hm xinT — e), Pn = pUn — 0) = lim p(nT — e), tn = nT (6) 

£-s>0+ e->-0+ 

is connected with the solution on the right-hand side of the kick a;(t„ + 0), p(t„ + 0) by Eq. (jl]), 
and the condition of continuity x{tn + 0) = x{tn — 0). The integration of (jl]) over the interval 

{tn - e, tn + e) gives 

p{tn + 0)=pitn-0)-KTG(Xn). 

Using notation ([6]), and the solution ([5]), we can derive the iteration equations 

Pn+1 =Pn- KTG{Xn). (7) 

Equations ([7j) are called the universal map. 

If G{x) = —X, then equations (^^ give the Anosov-type system 

Xn+l = Xn+ Pn+lT, Pn+1 = KTXn + Pn- (8) 

For G{x) = sin(x), equations ([7]) are 

Xn+i = Xn+ Pn+iT, Pn+1 = Pn - KT sin(a;„). (9) 

This map is known as the standard or Chirikov- Taylor map [1]. 



3 Fractional generalization of the universal map for 1 < 

In this section, a fractional generalization of the differential equation ([3]) is suggested. The 
discrete map that corresponds to the fractional equation of order 1 < a < 2 is derived. This 
map can be considered as a generalization of the universal map for the case 1 < a < 2. 
Let us consider a fractional generalization of ([3]) in the form 

CO 

^D'^x + KG{x)Y,^i^-f-n) =0, (1< a < 2), (10) 

where o-Df is the Riemann-Liouville fractional derivative [HI Ej [7] , which is defined by 

oDU = Dl o/r^x = — I--^ f ^^^^^^, (1< a < 2). (11) 

Here we use the notation D"^ = cP/dt"^, and qI^ is a fractional integration [HI El [7]. 

Proposition 1. The fractional differential equation of the kicked system / flOj) is equivalent 
to the discrete map 

Pn+l =Pn- KT G{Xn), (12) 

rpa—l " 

a;„+i = -— - Vpfc+iV;(n-A; + l), (1< a < 2), (13) 

1 [a] ^-^ 

where the function Va{z) is defined by 

V^{z) = z'^-' - (z - ir-' iz>l). (14) 

Proof. Let us define the variable ^(t) such that 

^D',--^ = x{t), (15) 

where qD^~°' is the Caputo fractional derivative 

fA-5^o/r'A'«^^/j^^, (0<2-„<l). (16) 



Using Lemma 2.22 of ^, we get 

oit"^Dt^^ = at)-m- (17) 

Then 

oD^x = D^ o/r-x = A' olt" ?A'-"e = D'Mt) - m) = A'e (18) 

Substitution of (HH]) and ([15]) into Eq. ([10]) gives 

n=0 

This fractional equation can be presented in the Hamiltonian form 



oo 

A^e + ^G(^Dr"OE'^(^-^)=0' (1<«<2). (19) 



oo 

^ + irG(o^A'-"OEKr-^) =0- (20) 

ra=0 

Between any two kicks there is a free motion 

rj = const, C, = Tjt + const. (21) 

The integration of (120]) over (t„ + e, tn+i — £) gives 

^{t„,+i-0)=^itn + 0)+Vn+iT, vitn+i-0) = r]{tn + 0). (22) 

The solution of the left side of the ra-th kick 

C.n = ^{tn - 0) = Urn ^{nT - e) , r]n = r]{tn - 0) = limr]{nT - e), t^ = nT, (23) 

£-S>0 £-5-0 

is connected with the solution on the right-hand side of the kick C,{tn + 0), ri{tn + 0) by equations 
([20]) . and the continuity condition 

^(tn + 0)=e(t„-0). (24) 

The integration of ( 120]) over the interval (t„ — e, t„ + e) gives 

7^(t„ + 0) = vitn - 0) - irTG(^A2,7"0- (25) 



Using notation (!23l) . and the solution (l2Ti) . we get 



r]{tn + 0) = r](tn+i - 0) = r7„+i. 
Substituting of 1^^ and (p5|) into (!22|) . we derive the iteration equations 

Vn+i = Vn-KTGi^Dl'^0- 
To derive a map, we should express the fractional derivative 

1 /•*" dr d^(r) 






through the variables (1231) . Using ^ = r^, we have 






r(«-i)7o (tn-r) 



'" r]{T)dT 

2-a' 



It can be presented as 



ra-l 



Cr,2-a. 



1 ^^ r'=+^ r7(r)cir 



where t^+i = t^ + 7" = (A; + 1)T, t^ = kT, and to = 0. 
For the interval (tfc,tfc+i), equations (|2T]) and (l23l) give 

?7(r) = r7(tfc + 0) = r/(tfc+i - 0) = %+i, r G (4, tfe+i). 



(26) 

(27) 
(28) 

(29) 



(30) 



(31) 



(32) 



Then 



*'=+! ri{T)dT 






2-a 



-fc+i 



Vk+1 



a — 1 



Vk+i / (tn - r)" dr 

'ifc 



^° ^2; = r]k^ 



z^ 



1- 



tn—tk + 1 



Vk+1 



tn—tk 
t-n—tk + l 



a — I 

[tn — tk) — [tn — tk+1) 



Vk+1 



T 



a-l 



a — 1 



{n - kY~' ~{n-k-l) 



Q-l 



(33) 



Using (a — l)T{a — 1) = r(Q;) and Eq. (!33l) . the fractional derivative (pTl) can be presented as 

o''Al-"e=^E^'^+i^"(^-^)' (1<«<2), (34) 

^ ^ k=0 

where 



14(2) = ^°-^ - (2 - 1)"-^ (z > 1). (35) 
As a result, Eq. ( l28l) takes the form 

^n+l = ^n + ?7n+lT, (36) 

W+i = r7„-irTG(— — 5^77fe+iV;(n-A;)j, (1< a < 2). (37) 



These equations define the fractional generalization of the universal map for variables {C,n,Vn)- 
The fractional equation flTOj) in the Hamiltonian form can be presented as 



°° t 



p + KGix)Y,s{^-n)=0, (38) 



n=0 



oDr'x = p, (l<a<2), (39) 

where we use qD" = D\ oD^~^. 
Equations flTSl) and fl34|) give 

^ ^ k=0 

Using equations fITTj) and fl39l) . we obtain 

p = oDr'x = Di o/r-x = Di o/r- ^A'^^e = A'(ew - ^(o)) = A'e (41) 

The definition of rj in fl20|) and Eq. (HTj) give 

P = i = V, Pn = Vn- 
As a result, equations (1H7I) . and (HU]) give 

Pn+l=Pn-KTG{Xn), (42) 



T 



a-l 



•^n+1 



^ ' A;=0 



Y,Vk+iVo,{n-k + l), (l<a<2), (43) 



where Va{z) is defined in (135!) . 
This ends of the proof. D 

Equations (H2!) and (H3l) define the fractional universal map in the phase space (x„,p„). 
These equations are generalization of the map (171). 

Note that the form of Eqs. fH3l) is defined by both equations fl38|) and fl39|) . Equations P3|) 
cannot be considered as an iteration representation of Eq. ( 139|1 only. If we use the other form 
of Eq. ([38]), then Eq. (USD is changed. 

Let us consider some examples of fractional universal map f l36|) . fl3l 



Example 1. Let us prove that the fractional universal map for a = 2 gives the usual 
universal map. Substitution of Eq. (l36l) in the form 



r7fc+i = -(a+i-6) (44) 



into the iteration equation (P7|) gives 



r/„+i = T^n - KTg(— ^ ^(e;.+i - 6)V;(n - A;)). (45) 

^ '^ fc=0 

Then, the fractional map fl36|) . f l37|) is defined by 

Sra+l = Sn + Vn+lT, 

r7„+i = r/„ - iTTG^-— ^^(^^+1 - ^k)V^{n - fc)), (K a < 2). (46) 



For a = 2, we have Va{z) = 1, and 



n — 1 



J2(^k+i-^k)=^n-^o. (47) 



fc=0 

Then equations fj45|) gives 



7]n+l=Vn-KTG(^^n-^o)- (4J 



9 



As a result, Eqs. ( H6l) with a = 2 give the usual universal map l^^ for the variables Xn = ^n — ^o, 
and pn = r]n. 

Example 2. If G{x) = -x, then Eqs. ll36l). fl37|) are 



^n+l — 'Cn. + Vn+lT, 



T 



n-l 



1171+1 = 1171 + K——-^r]k+iVa{n-k), (1 < a < 2), 



(49) 



k=0 



where Va{z) is defined in (l35l) . For the space (x„,p„), this generalization can be presented as 



rpa 



n-l 



Pn+1 =Pn + K——y^^pk+iVa{n - A; + 1) 



fc=0 



T 



a-l 



X 



n+1 



T(a 



^Pk+iVa{n - k + 1] 



(50) 
(51) 



k=0 



Equations flHI]) and 05 ip define a fractional generalization of the Anosov-type system. 



Example 3. If G{x) = sin(x), equations (j36|) . (I37j) give 



^n+l — 'Cn + Vn+lT, 



Vr 



rpa-1 "-1 , 

+1 = ?7„ -/CT sinf— — ^r7fc+iVa(ra- A;)j, (1 < a < 2). 



(52) 



This map can be considered as a fractional generalization of the standard map. The other 
possible form of Eq. fl52|) is 

Sn+l = 'Cn + ^n+l5 
1 



^n+1 = Vn-K sin(— - ^(^fc - ^fe_i) (n - A;)" ^ - {n - k - 1) 



k=l 



\a-l 



where we use Eqs. flUl) . (1551) . T = 1, and 1 < a < 2. For (x„,p„), equations 
G{x) = sin(x) give 

Pn+i = Pn - KT sin(x„), 



(53) 
with 

(54) 
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T 



a-l 



•^n+1 



^ ' fe=0 



Y,Vk+iVo.{n-k + l), (l<a<2). (55) 



These equations define a fractional standard map on tlie pliase space, wliicli can be called the 
fractional Chirikov- Taylor map. 

Example 4. The fractional generalization of the dissipative standard map [HI [19] can be 
defined by 

in+l = in + Vn+l, (56) 

r]n+i = bVn-K sin(— ^ J^^^^ - ik-i) [{n - k^-' - {n - k - 1)"-^]). (57) 

^ ' fc=i 

For 6 = 1, we get the fractional standard map (1531) . 



4 Fractional universal map for a > 2 

In this section, a fractional differential equation f lTU]) is used for a > 2. The discrete maps 
that correspond to the fractional equations are derived. These maps can be considered as a 
generalization of the universal map for the case a > 2, i.e., the fractional universal map fl36|) . 
( 137|) can be generalized from l<a<2toa>2. 
Let us consider the fractional equation 

oo 

oD'^x + KG{x)^5{--tt\ =Q, (m-l<a<m), (58) 

n=0 

where o-D" is the Riemann-Liouville fractional derivative of order a, m — 1 < a < m, which is 
defied [S [6l [7] by 

oDtx = DT ^ir^^x = —^ " / ^ ^^^^"^ ^ , (m-l<a<m). (59) 

Here we use the notation D^ = d'^/dt"^, and o-^t"~" is a fractional integration [HI [Gj [7|. 



11 



Proposition 2. The fractional differential equation of the kicked system ( fjgj) is equivalent 
to the discrete map 



n—l »n— 3 



Xr, 



ly E(E ii Pn'V^'^i^ -k) + Z^^r^P^+i' V-'-'{n - k)), (60) 



r(a-m + l)^V^ n -" - ^ ^ {m-2)V 

m—s—2 rjii r-prn—s—l 

Pn+l=fn+ E 7rPn^'+ (^„,_i), P:r+/. (^ = 1, .-, ^ - 2), (61) 

Pn+i' = Pr' - ^TG(x„), (m - 1< « < m), (62) 

where the functions V^'\z) are defined by 

Kr''W= r (^-y)y-™rf2/= / y'iz-yT^^'dy, / = l,...,m-2, (63) 

Jz-l JO 

and l</<'m — l<a<'m, where l,m gN. 

Proof. Let us consider equation (fTOl) with m — 1 < a < m. Using ^ = ^(t), such that 

^D["-"^ = x(t), (m-l<a<m), (64) 

we obtain 

oA"^ = oA" o^^A^^-^e = A" o/r"° ^A'"""^ = A"e, (O < m - a < l). 
Equation (ITOl) can be presented as 

oo 

Dre + i^G(o^A'"""0 E'^lr - ^) = 0' (^ - 1< « < ^)- (65) 

n=0 

Let us define 

r]'{t) = Dtm, {s = 0,...,m-l). (66) 

Then the Hamiltonian form of Eq. (!65|) is 



oo 



r]^ = n'+\ (s = 0,...,m-2 

n=0 

12 



r-' +KTG{ ^Dr^v') E '^ (r - ^) = °- ^^^^ 



For t E {tn + 0, tn+i — 0), equation ( l65l) is 

Dn = 

and the solution can be presented as 

m— 1 

eW = 5ZQ(t-t„)\ (m>3) (68) 

z=o 

Substitution of fl68|) into fl66|) gives 

m— 1 

r7^(t) = ^ Q/(/ - !)...(/ -s + l){t- tJ-\ (69) 

l=s 

For t = tn, v/e have 

r/^(t„ + 0) = C,s!, 

Then 

a = ^r/'(t„ + 0) (70) 

Using Eq. f irUjl . the relations 

r^^(t„ + 0) = 7]'itn - 0) = r/:, (s = 0, ..., m - 2), 

and 

^-i(t„ + 0) = r^™-i(t„+i - 0) = C+-i\ 

we present equation (1691) in the form 

m-2 



1 1 

^'W = E7^ ^i(i-in)'"^ + 7 7T;Vn^iit -tnr-'-\ is = 0,...,m-2). (71) 

■^ (/ — sj! (m, — s — Ij! 

/(/-!)...(/ -s + 1) 1 



Equation (ITT]) can be rewritten as 



^^W= E 7r^:^'(^-^-)' + 7 ^^TTC+i'(t-U"^-^-S (s = 0,...,m-2). (73) 

■^^ i! m — s — 1 ! 



i=0 

13 



For t = tn+i, this equation gives 

m— s— 2 



'n+i= E Wn^'^' + l ^^TtC+-i'T— \ (. = 0,...,m-2). (74) 



As a result, the iteration equations are 

m—s—2 



Vn 



1= E k^'^'+ ,^_! ni ^""^+i'^""^''' (. = 0,...,m-2), (75) 



/! '" (m-s-iy. '"+ 



Vn+l = Vn-' - KT G{ ^D^rv')- (76) 

The Caputo fractional derivative 

can be presented as 






where t^ = fc^- 

For t G (tfc^^fc+i), equation fITH]) is 

m— 3 -. ^ 

/=o ■ I J- 

Substituting of fl75]) into fl77|) . and using 

.. (i„ -*)'""" it (n-r)"-" 

we obtain 

ftfc+i Til/'fWf "1-3 rpa+l-m+1 rpa~l 

where 






V"'Hn-k)= -^ -^dz, l = l,...,m-2. (80) 
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These functions can be defined by 



v"'Hz) 



z-1 



{z-y)Y-^dy= y'{z-yr-^dy, / = !,... ,m-2, 
Jo 



wfiere 1 < I < m — 1 < a < m, 1,171 &N. 
Let us use equation flM|) in the form 



X. 



C j-^m—a 



n — 



Ar% Pn = v:, (s = i,...,m-i). 



As a result, we obtain 



n-l m.-3 rpa+l-m+1 



Xr, 



T(a-m + l) ^^\^^ l\ 

^ > k=0 1=0 



T 



a-l 



E ( E ^^^ Pn^'y^'in -k) + J^^;^^Pn^i V^^--\n -k)), (82) 



m-s-2 rpl 



Pn+1 ~ Pi 



+ E Tp;"' 



T 



rn—s—l 



1=1 



{m — s — 1)! 



p^~l, (s = l,...,m-2), 



p;r+/ = Pr' - ^TG(x„), (m - K « < m). 



(83) 
(84) 



This ends of the proof. D 



Remark 1. Equations fl82l) - fl8^ define a fractional generalization of the universal map for 
a > 2. For G{xn) = Xn, equations (I82 l) -( l84l) define fractional Anosov-type system with a > 2. 
For G{xn) = sinx^, we have the fractional standard map for a > 2. 



Remark 2. The functions that are defined by integrals ( IHTI) can be expressed in elementary 
functions. For example, equation (l80l) with / = gives 

1 



V7'°(n - k) 



a — m + 1 
This equation for m = 2 can be presented as 



{n — k) 



a—m+l 



{n-k-l) 



a—m+l 



5) 



1 



For / = 1, Eq. ([80]) gives 



V^^\n-k) 



(a — 771 +!)(« — m, + 2) 



V^^\n-k) = -V^{n~k). 

a — i 

{n-kY-'^+'^-{n-k-lY-'''+^{n-k+a-m+l) . (86) 
15 



Remark 3. Note that the function ( IHTl) can be presented through the hypergeometric 
function (section 2.1.3. of [21j) F{a,b,c; z) by the relation 

Kr'(^) = ii + X)z^~. n^ - a, / + 1, / + 2; z-'l (87) 



where 



and we use 



r(/ + i) i\ 1 



V{l + 2) (/ + !)! / + !' 

5 Fractional kicked damped rotator (FKDR) map 

In this section, a fractional generalization of the differential equation for a kicked damped 
rotator is suggested. The discrete map that corresponds to the fractional differential equation 
is derived. 

Let us consider a kicked damped rotator [2]. The equation of motion for this rotator is 



x + qx = KG{x)^5(t-nT\ (88) 

n=0 

It is well known that this equation gives [2] the two-dimensional map 

yn+i = e-'''^[yn + KG{xn)l (89) 

1 - e-«^ 

Xn+l = Xn^ [yn + KG{Xn)\. (90) 



Let us consider the fractional generalization of equation (18811 in the form 

oo 

oD^x-qoD^^x = KG{x)^5{t-nT), (91) 



n=0 

where 

gGM, l<a<2, /3 = a-l, 
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and o-Df is the Riemann-Liouville fractional derivative |5l El [7] defined by ( TTTl) . Note that we 
use the minus on the left-hand side of Eq. fl9T]) . where q can have positive and negative values. 



Proposition 3. The fractional differential equation of the kicked system ( fglj) is equivalent 
to the discrete map 

Pn+l = e'''^ \pn + KG{Xn) 
rpa—1 " 



Xn+1 



r(a-l) 



^Pk+iWaiqT,k -n-1] 



k=0 



where the functions Wa{a,b) are defined by 



Wa{a, h) = a^-"e'^(^+^) [r(a - 1, ab) - T{a - 1, a{b + 1)) 



andT{a,b) is the incomplete Gamma function: 



T{a,b)= / y'^-'e-ydy. 

Jb 

Proof. Let us define the variable ^(t) such that 



^D',--^ = x{t), 



where qD^ " is the Caputo fractional derivative flT6|) . Using 



we obtain 



^j2-a Cjj2^a^ ^ ^(^) _ ^(q)^ (Q < 2 - a < 1), 



oD^x = Dl ol't-'-x = Dl o/r- '^Dj-^i = D^im - m) = D^, 



and 



Substitution of ([98]), ([99D and ([96]) into Eq. ([H]) gives 

oo 



n=0 
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(92) 
(93) 



(94) 



(95) 



(96) 



(97) 



(98) 



oD^x = D] oll-'x = D] olf^x = D] ,P,-- ^^^-^^ = Dl{m - ^(0)) = D]i. (99) 



(100) 



The fractional equation fllOOp can be presented in the Hamiltonian form 

oo 

Between any two kicks 



n=0 



■f] — qrj = 0. 
For t G (t„ + 0, t„+i — 0), the solution of equation (11021) is 

r]{tn+i - 0) = r]{tn + 0)e^^. 

Let us use the notations 

Cn = x{tn-0) = \im ^{nT - e), rj^ = rj{tn - 0) = limrj{nT - e), tn = nT. 

For t G (t„, — £, tn+i — £), the general solution of fllOip is 

oo „t 

. n Jt„—e 

Then 

r]n+i = e" 

Using (11061) . the integration of the first equation of (IIOII) gives 

l-e"'^ 



m=0 



olT 



Vn + KG{l,Dt;-0 



Sn+1 — Cn 



r^. + KG(o^A„ 



:ior 



(102) 



(103) 



(104) 



(105) 



(106) 



(107) 



Let us consider the Caputo fractional derivative from Eqs. (I106p . (I107p . which is defined by 



C Tt2-a 



\^Dtr^=oir'DU 



Using C, = T], we have 



dr d^{ 



r 



T{a-l)Jo (t„-r)2-" dr 



(0 <2-a< 1). 



(108) 



° *" 4-oi, ''r(a-l)7o (t„-r)-i- 



*" rj{T)dT 



(109) 
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It can be presented as 






^ n—l 



■'=+^ rj{T)dT 



r("-l)t^4 (^n-rp-- 



(110) 



where t^+i = t^ + ^ = (^ + l)^? ^fc = ^7", and to = 0. 
For r G (tfc,tfc+i), equations f llOSp and fll04p give 



Then 



'- „(r).. .„^j"%,.-.«,(,„ _,)<.-.,, 



(t„- 



T 



tk 



"tn—tk 
'in— tfe + l 



in ifc 



z"-"e-''"t/2 



tn—tk + 1 



r/.+ig^-"e«("-'^-^)^ 



g(in-tfc) 



y^-^e-ydy. 



(112) 



<j(tn— tfc+l) 



As a result, Eq. (I112p gives 



'fc+ 



1 ^(t)c^^ _ l_a „(„_fc_i)T 

(t„ - r)"- - ^'^^^^ ^ 



r(a - 1, q{tn - tfe+i)) - T{a - 1, g(t„ - 4)) . (113) 



Here r{a,b) is the incomplete Gamma function [20] : 

/oo 
y^-^e-^dy, a,6G 



which can be defined by 



r(a, b) = r(a) iFi(l, 1 + a; -6), 

a 



where iFi is the confluent hypergeometric Kummer function 

(a)fc2;^ 



iFi{a,c;z) = ^ 
Here (a)^ is the Pochhammer symbol 



k=0 



ic)k k\ 



(114) 



(115) 



(116) 



(a)fc = a(a + l)...(a + fc- 1), A; G N. 



(117) 
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Using fITTnll and flTT3ll . we get 



?Al'"e = ^J^^. Y.'^k+iWMT.k-n), (1< a < 2), 



fc=0 



where 



Wo,{a, h) = ai-°e"(''+^) r(a - 1, afe) - r(a - 1, a(6 + 1)) 
Substitution of (ITTSD into (IT06|) and (ITOTl) gives 



Vn+i = e 



qT 



r(^ - 1) 



n-l 



Sn+1 — Sn 






fc=0 



n-l 



^" + ^^(-^^^ 5^r7fc+iW^«(gT, fc 



r(a-l) 



n 



(118) 



(119) 



(120) 



:i2ii 



fc=0 



Equations (11201) and (11211) can be rewritten as 



n-l 



Vn+i = e 



qT 



Vn + KG{ ^^I^^^J2vk+iWa{qT, k 



r(a-l) 



n 



k=0 



, _ , 1 - e-"^ 

^n+l — C,n H Vn+1- 

q 



:i22) 



(123) 



These equations can be considered as a fractional generalization of kicked damped rotator map 

for {^n,Vn)- 

Equation f l9T|) can be presented in the Hamihonian form 



oD^-^x = p, p-q oD^x = KG{x) ^ 5(t - riT) . 
For {xn,Pn), equations (!96|) . ( |98|) . (!99|) give 



(124) 



n=0 



Pn ^n; 2Jj; 



?A^-"e 



As a result, we have 



p„+i = e"'^ ^pn + KG{Xn) 



T 



•a—l 



Xn+l = :^-———'^Pk+lWa{qT, k - 71 - 1) 



r(a-l) 



(125) 
(126) 



fc=0 
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where Wa{a, h) is defined in f lllQp . This ends of the proof. D 

These iteration equations fll25p and (11261) define the fractional kicked damped rotator map 
( l88ll in phase space (x„,p„). 



6 Conclusion 

There are a number of distinct areas of physics where the basic problems can be reduced 
to the study of simple symplectic maps. In particular the special case of two-dimensional 
symplectic maps has been extensively studied. Under a wide range of circumstances such maps 
give rise to chaotic behavior. The suggested fractional maps can be considered as a fractional 
generalization of the symplectic map that is derived from kicked fractional differential equations. 
We can suppose that fractional discrete maps can be connected with some generalization of the 
symplectic structure. Fractional generalization of symplectic structure [22] can be defined by 
using fractional differential forms. 
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